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Abstract 

We consider the compressible Kraichnan model of turbulent advection with small 
molecular diffusivity and velocity held regularized at short scales to mimic the ef- 
fects of viscosity. As noted in 0, removing those two regularizations in two opposite 
orders for intermediate values of compressibility gives Lagrangian flows with quite 
difi'erent properties. Removing the viscous regularization before diffusivity leads to 
the explosive separation of trajectories of fluid particles whereas turning the regu- 
larizations off in the opposite order results in coalescence of Lagrangian trajectories. 
In the present paper we re-examine the situation first addressed in 6 in which the 
Prandtl number is varied when the regularizations are removed. We show that an 
appropriate fine-tuning leads to a sticky behavior of trajectories which hit each other 
on and off spending a positive amount of time together. We examine the effect of 
such a trajectory behavior on the passive transport showing that it induces anoma- 
lous scaling of the stationary 2-point structure function of an advected tracer and 
influences the rate of condensation of tracer energy in the zero wavenumber mode. 
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1 Introduction 

Certain qualitative features of transport of scalar or vector quantities by highly turbulent 
flows may be understood using drastically simplified models where the back-reaction of 
the transported quantity on the velocity dynamics is ignored (passive advection) and the 
hydrodynamical velocities are replaced by a simple synthetic random velocity ensemble. 
In recent years much attention has been attracted by the study of passive advection 
in a Gaussian ensemble of time-decorrelated velocities, first considered by Kraichnan 
[T5| and Kazantzev |14,. The model has led to a new understanding of the statistical 
intermittency of advected quantities and of the interplay between fluid compressibility 
and transport properties, see FT and references therein. The passive advection by a d- 
dimensional velocity field z/(t, r) is intimately related to the behavior of the Lagrangian 
trajectories (i.e. trajectories of fiuid particles) perturbed by noise. Such trajectories 
satisfy the stochastic ordinary differential equation 

^^R{t)^v{t,Rit)) + V2^rf{t), (1) 

where f]{t) is a d-dimensional white noise independent of velocity (and of noises of other 
particles) and k > is the (molecular) diffusivity. The advection of a passive tracer 
0{t,f) is governed by the advection-diffusion equation 

dtO + V- V0 - KV'^e ^ f , (2) 

where f{t,r) is a source or forcing term. In the absence of forcing (i.e. for / = 0) the 
passive tracer is carried along the particle trajectories so that 

e{t,r)^e{Q,R{Q-t,T)). (3) 

Here R{t' ] t, r) denotes a particle trajectory that passes through point r at time t and the 
over line stands for the average with respect to the noise fj{t). Note that the forward tracer 
evolution is described by Lagrangian trajectories going backwards in time. Similarly, in 
the presence of the source, the tracer evolves according to the equation 

e{t,r)^e{{),R{0-t,r))+ f f{s,R{s;t,r))ds. (4) 

Jo 

This way, in random velocities, the statistics of advected quantities like 9 is linked to the 
statistical properties of (noisy) trajectories. The present article exploits this relation in 
yet another situation in the Kraichnan model of turbulent advection. 
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1.1 The Kraichnan model 



The Kraichnan model describes advection by a stochastic velocity field v{t, r) with Gaus- 
sian mean-zero statistics, stationary and decorrelated in time, homogeneous and isotropic 
in space and with fixed compressibility degree p characterizing the relative strength of 
the incompressible and potential components of the velocity. As in any centered Gaus- 
sian ensemble, the statistics of velocities is totally characterized by the 2-point function 
(^(t, r) ® v{t' , r ')) . In what follows, we shall need the statistics of the equal-time velocity 
differences that is determined by the reduced expression 

EE Sit-t')d{r-f'). 

The rank 2 tensor d{r' — is the second order spatial structure function of the velocity 
field. The rank 2 tensor P(fc, p), invariant under rotations and of trace 1, is taken to be 



with < p < 1. The coefficient p is equal to g'g'j^'^ (summation convention!) for all 

r 7^ 0, so that it is meaningful to characterize p as the ratio l^g'^'jaj even though the 
velocity field is not smooth for ^ < 2. The value p = corresponds to an incompressible 
velocity field, whereas p = 1 to a potential (e.g. irrotational) one. Note however that 
if the physical space is one dimensional (case d — \) then necessarily p = 1 because a 
one-dimensional velocity field is always potential. One simply takes P = 1 in this case. 

The spatial velocity structure function d{r) scales in |r| with power ^. We shall take ^ 
between and 2. This choice guarantees that the typical velocity realizations are non- 
Lipschitz. More exactly, they are Holder-continuous with exponent i^' for any ^' < ^ 
(the realistic turbulent velocities are believed to have Holder exponent < in the limit 
of infinite Reynolds numbers JH!)- Note that the (positive) constant Dq in ^ has the 

J. . (length)^"? 

dimension — . 

time 

1.2 Regularizing effects 

Real flows are always regularized at small scales by viscous effects. That is, the velocity 
field is smooth and the power law scaling with ^ < 2 is observed only at distances much 
larger than the viscous cutofT scale li, that becomes very small only for high Reynolds 
numbers (if the integral scale of turbulence is kept fixed). The small r behavior cx \r\^ of 

d{f) above comes from the the slow decay at large k of its Fourier transform (the term 
under the integral sign in Q). We shall mimic the effect of viscosity in the present model 

by introducing an ultraviolet regulator and replacing d{r) by 

where the function / is taken positive, smooth, decreasing, fast decaying at infinity and 
with /(O) = 1, /'(O) = 0. The viscosity v itself may be defined as 

u - D^li (7) 

which is the only combination of Dq and ly of dimension ^'^".^^-^ . 

Another regularizing effect in real flows comes from finite molecular diffusivity. The 
tracer field 0{t^ r) passively transported by the flow will be smoothed by diffusion, see 
0. In typical velocities, this smoothing effect becomes important at scales much smaller 
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than the diffusive cutoff scale Z^. Beneath this scale diffusion dominates advection. We 
modeled molecular diffusivity by adding white noise terms to the fluid particle velocities 
(independent ones for different particles, see (0). In the Kraichnan model, the scale Z^ 
beneath which diffusion dominates advection may be expressed in terms of k, Dq and l^,, 
with = corresponding to k = 0. The expression may be inverted to calculate k in 
terms of Dq,Ii^,Ii,. Later on we shall specify such a relation in the case when <C li,, 
see (|28f) . This will provide an expression for the Prandtl number Pr = ^. 

1.3 Statistics of fluid particles 

Even in a simple random velocity ensemble, the statistics of the Lagrangian flow may 
be quite complicated. It may be studied by looking at the joint A'^-particle probability 
density functions (PDFs) defined by 

N 

Here, as before, the overline denotes the average over the (independent) white noises 
and ( ■ • • ) stands for the velocity ensemble average. In the Kraichnan model, due to 
the temporal decorrelation of velocities, the PDFs (O are Markovian and they define a 
consistent hierarchy of A''-particle stationary Markov processes that contains the complete 
information about the statistics of the Lagrangian flow and of the velocities themselves 
|17| . In this paper we shall be interested uniquely in the behavior of the separation of a 
pair of Lagrangian particles. The main object of our interest will be the PDF of finding 
their time t separation equal to r, given that their time t' separation is equal to tq: 

P*'*'(f,f') = (^6{f' ~ Rit';t,f) + R{t';t,6))y 

In stationary, time-reversal invariant velocity ensembles, such as the Kraichnan model, 
P*'* depends only on |i — t'| and we shall use the notation pl*~* I. 

The PDF P*'* (r, r") governs the free decay of the 2-point function of the passive 
tracer, evolving according to equation with / = 0. If at time zero the tracer is 
distributed independently of the velocity field and the trajectory noises, with a homoge- 
neous 2-point function {9{0,r) 6{0,f')) — P(0,r' — r), then at a later time t its 2-point 
function is given by 

P(t, f)^ J P(0, f') P''°{r, r') dr' (9) 

as follows from Similarly, if the scalar source / is a random field with mean zero 
and 2-point function 

{f{t,r)f{t',r'))=5{t-t')x{r~r') (10) 

and it is independent of the velocity field, the trajectory noises and the initial scalar 
distribution, then the evolution of the tracer 2-point function is described by the relation 

F(t,r)^ j F{Q,r')P^ °{r,r')dr' + ds J xir') P''" {r,r') dr' (11) 

as follows from Q by taking averages. 

Our aim is to describe the particle pair separation in Kraichnan velocities at large 
scales, much larger than the cutoff scales li, and Z^- The effective description will depend 
on the Prandtl number and on the scales involved. Instead of maintaining the cutoffs 
finite, we want to give an effective large scale theory without cutoffs but with some 
specific boundary condition at vanishing separation. Such a theory should give rise to 
the same behavior at the large scales of interest. An alternative way of doing this is by 
taking the cutoffs li, and 1^ to zero in a fashion that preserves large scale behavior and 
to examine the boundary condition that arises in this limiting process. 

^In the Kraichnan velocity that is white in time, Eq. should be interpreted as a Stratonovich 
stochastic differential equation, see e.g. Sect. II. C. 3 of 
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1.4 Lagrangian dispersion 

As first derived in 15 , the time evolution of the scalar 2-point function of the tracer in 
the Kraichnan ensemble of velocities is governed by the differential equation 

dtF = M.^^F 

where M^^k, is a partial differential operator which can be written in terms of the velocity 
structure function d and diffusivity k as 

Mu^f, = dij{r; l„)didj + 2kV^. 

It follows that the separation PDF takes the heat kernel form 

P*(fo,f) = e-^-(ro,r). (12) 

As we shall be only interested in distances between two particles, usually called the 
Lagrangian dispersion, and not in the angular distribution of the particle separation, 
we may project H12|l to the rotationally invariant sector. This is done by restricting 
the action of operator A4u,k to functions of the radial variable r = |r| only. A4u,k is a 
rotationally invariant operator and it maps functions of r to functions of r. We shall 
denote by M^^^ its restriction to such functions. It is easy to show that 

M,,^ = Do[fi{r)d^,+f2{r)dr] (13) 

with the coefficient functions fi^2(j) = fi,2(j']ly,lK) expressed in terms of the spatial 
velocity structure function d by 

Dofiir)^d,j{f;l,f-^ + 2K (14) 

Dof2{r) = id,,(f; mS,, - ^) + (15) 

for any r such that |r| = r. For detailed calculation of the values of /i and /2 see 
Appendix E It turns out that, for r much larger than the cut-off scales 1^,1^^ the 
function /i is asymptotically proportional to and the function /2 to r^~^ . We choose 
the constant Dq proportional to I?Oj see (jSJ, so that r^^ fi{r) — > 1 when r — > oo. 

In order to turn M^^^ into a true operator on the half-line M+, we also have to specify 
its domain ©(M^ ^)- This boils down to the choice of the boundary condition to impose 
at zero. The domain of A^i/^k, as long as k is positive, is some locally Sobolev space. 
For its rotationally invariant elements this implies that their gradient vanishes at r = 
resulting in the Neumann boundary condition dr\r=o = for M^^^. This way equation 
()12(l reduces in the rotationally invariant sector to the heat kernel expression for the 
particle distance PDF 

P*(ro,r)=e**^-''(ro,r), (16) 

where Af^^„ is the second order differential operator on the half-line of r > given by 
(|13(l . taken with the Neumann boundary condition at r = 0. We may view the right hand 
side of (11611 as the transition probability density of a Markov process r{t) that describes 
the distance between two noisy Lagrangian particles in the random flow and treat the 
underlying diffusion process with the adapted tools, see Appendix ^ Alternatively, we 
may interpret H16(l as describing kernels of a semigroup of operators and apply to it the 
usual analysis toolkit. Both methods will be developed. 

1.5 Initial classification of boundary behaviors 

The possible boundary behaviors at the origin of the one-dimensional diffusion process 
r{t) describing the inter-particle distance in the absence of regularizations (i.e. for 1^,1,^ = 
0) may be easily classified j6j. In this case the generator (fl^ of the process reduces to 
the operator 

M = Mo,o = Dor^ [dl + ^dr) (17) 
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with 

= (18) 
see Appendix E Let us note in passing for later use that we have the imphcation 

< 1 =^ o^a^p < 1 ■ 

For d>2 this is a simple consequence of H18|) . In the case d = 1 recall that p = 1 and 
then a^^i = for any value of ^. Below, whenever the value of the compressibility degree 
p is understood, we shall write instead of the more unwieldy a^^p and sometimes 
simply a if no confusion may arise. 

In the case ^ < 2 corresponding to spatially rough flows, we can introduce the new 
coordinate 

u = u(r) = ^^(7'/V-3. (19) 
In this coordinate the generator H17|l becomes 

+ (20, 

and it describes the Bessel process (see [2] for a brief resume) of parameter —6^ or 
effective dimension doff — 2(1 — 6j) where 

'. = ^- ™ 

If deff is a positive integer then the corresponding Bessel process describes the behavior 
of the norm |VF(<)| of the deff-dimensional Brownian motion W{t). 

The general theory of the boundary behaviors of one-dimensional diffusion processes 
was laid down by Feller in [Hj. For the Bessel process, the boundary behavior at zero is 
well known pil2U| and it depends on the parameter or equivalently the effective dimension 
of the process. For dog < 0, zero is an exit boundary point (the realizations may arrive 
at zero in a finite time, but they cannot come back from zero nor start there). For 
< deff < 2, zero is a regular boundary point (the realizations can go to and leave zero 
in finite time). For deff > 2, zero is an entrance boundary (the realizations may start from 
zero, but no realization goes to zero in finite time). The different boundary behaviors 
enumerated above correspond respectively to the regimes in the Kraichnan model that 
were coined strongly compressible (p > ^), of intermediate compressibility 

(^ + i < p < ^) and weakly compressible (p < ^ + \), see [TTIIK]. 

In the limiting case ^ = 2 corresponding to spatially smooth fiows, we may introduce 
the new coordinate u = Inr. In this coordinate the generator H17() becomes 

dl + (a2 - l)du 

and it describes the one-dimensional Brownian motion viewed in the frame moving with 
speed — (a2 — 1) (i.e. the Brownian motion with a drift). In this case r = (corresponding 
to M = — oo) is a natural boundary point with no realizations that attain it or leave it 
in finite time. The quantity (02 — 1)Dq^^ — i+tp -^o^'^ Lyapunov exponent of the 

Lagrangian flow and it is non-negative in the weakly compressible regime p < j, and 
non-positive in the strongly compressible one p > j, vanishing at their common point 
p — J. We sum up this classiflcation in Tabled The first three lines correspond to rough 
velocity fields with £, < 2. The last two correspond to smooth ones with ^ — 2. Note 
that for d > 4 the weakly compressible regime extends to the whole interval < p < 1. 

The relationship between the boundary behavior of the one-dimensional diffusion 
process r{t) describing the distance between Lagrangian trajectories, on the one hand, 
and the different regimes of the Lagrangian flow, on the other hand, may be explained 
intuitively. If the realizations can go to zero but cannot come back, then they are trapped 
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flow 


conipressibihty 


characterization 


implications 


boundary 


rough 


weak 


P<^ + l 




entrance 


intermediate 


^ + 5 < P < F 


1 > a<; ,^ > f — 1 

1 > Go 


regular 


strong 


P > 4- 


f — 1 > fl;; 


exit 


smooth 


weak 


P < T 


02. p > 1 


natural 


strong 


P> 3 


02, p < 1 



Table 1: Regimes of the unregularized Kraichnan model 



there. This implies coalescence of Lagrangian trajectories that characterizes the strong 
compressibility regime, as first observed in . If realizations can go to and come back 
from zero, then one has to specify further the boundary behavior. There may be (and, 
as we find, there are indeed) different possible behaviors of Lagrangian trajectories when 
they meet. This is what happens for intermediate compressibility, as first noted in 
If realizations starting from inside the half-line can never reach zero then no trapping is 
possible. In this situation, corresponding to weak compressibility, the realizations may, 
however, enter from zero if ^ < 2 meaning that the Lagrangian trajectories separate 
explosively leading to a non-deterministic Lagrangian flow . For ^ = 2 the realizations 
can neither collapse to zero nor explode from it. In smooth flows, different trajectories 
remain separated indefinitely in the future and the past, although they may behave 
chaotically, with small separations growing exponentially in time. 

1.6 Main result 

In any case, the regulated generator M^^f^ of Eq. taken with Neumann boundary 
condition at zero gives rise to a diffusion process r{t) on the half-line describing the 
statistics of the Lagrangian dispersion. We would like to know how this process behaves 
in the limit when 1,^,1^ — > 0. 

The cases of the weakly and strongly compressible regimes where the limiting disper- 
sion process does not depend on the way in which the limit is taken seem to have been 
treated adequately In these cases only one process is possible in the limit and it 

corresponds either to instantaneous reflection of trajectories when they meet (for weak 
compressibility) or to coalescence of trajectories upon the first hit (for strong compress- 
ibility) . 

In ref. [Sj |S] it was argued that, in the intermediate regime, the possible limiting 
processes are described by the unregulated generators M of Eq. H17() with boundary 
conditions at zero depending on the way the limit 1^,1^ is taken. Our analysis 
confirms this picture but differs from 6 in specific conclusions (the nature of the limiting 
boundary conditions and the form of the limiting procedures that lead to them). Upon 
the fine-tuning H29|l of l^, and l^, we find in the limit the so called sticky or slowly 
reflecting boundary conditions parametrized by the amount of "glue" A e [0, oo]. 

The sticky behavior can be pictured intuitively in the following way. For simplicity 
let us consider just the case of Brownian motion on the half-line. Concerning ordinary 
reflected Brownian motion, it is a well known fact that, if it hits the end-point of the 
half-line at some instant t, then for any instant t' > t it will hit again between t and 
t' an inflnite number of times almost surely, but still the cumulated time spent at zero 
is zero almost surely. Slowly reflected Brownian motion returns also almost surely an 
infinite number of times after the first hit of the boundary, but the cumulated time spent 
at zero (technically speaking the Lebesgue measure of those instants in [t, t'] when the 
particle is at zero) will be proportional to the time spent in an infinitesimal neighborhood 
(0, dx) of (not containing 0!), the proportionality constant being X/dx. Remarkably, 
however, for any X < oo, the trajectory will almost surely never remain at for an 
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uninterrupted strictly positive time-interval. The two extreme values of A are A = and 
A = cxD. In the first case zero time is spent at zero and, in fact, one recovers the ordinary 
(instantaneous) reflection. The second case with infinite time spent at zero corresponds 
to absorption (more precisely, adsorption) at zero. Transposed to the case of Lagrangian 
dispersion, this means that particles upon hitting each other either have no interaction 
(case A = 0), or they will tend to stay together for some positive amount of cumulated 
time (case < A < 0), or they stick together permanently (case A = oo). 

We shall give strong general arguments in favor of the convergence of the particle 
dispersion to the sticky process and shall prove the convergence for specific quantities. 
To analyze the effect of stickiness on the passive transport, we shall construct and analyze 
the behavior of the transition probabilities of the sticky dispersion process that govern, 
in the appropriate 1^,1^ limit, the evolution of the 2-point correlation function of 
advected tracer in the homogeneous and isotropic situation. It appears that the presence 
of stickiness induces anomalous scaling of the stationary 2-point structure function of the 
forced tracer and that it influences the rate of tracer energy condensation in the constant 
mode. For the two extreme cases of instantaneous reflection and absorption, not only 
the particle dispersion process but the entire family of the consistent A^-particle Markov 
processes has been constructed in ^| El- In the sticky case, the existence and the 
uniqueness of such a construction is still an open problem and is not discussed further 
here. Let us only note that ref. |18j may be interpreted as providing such a construction 
for the limiting case of the one-dimensional flow with ^ = 0. 

Let us recapitulate the situation that we are dealing with. We are studying the inter- 
mediate compressibility regime of the Kraichnan model of passive advection with velocity 
field smoothed at scales smaller than l^, and with molecular diffusivity n dominant on 
scales smaller than l,^. We shall be looking at the 2-point function of the advected 
tracer in the isotropic sector or, alternatively, at the PDF of the distance between two 
noisy Lagrangian trajectories. We want to know how those objects evolve in the limit 
IvJk. 0. The remainder of the paper is organized as follows. In Sect. [3 we shall discuss 
the probabilistic aspects of the problem studying the one-dimensional diffusion process 
describing the inter-particle distance. The analytic approach based on the spectral anal- 
ysis of the generator of the process will be developed in Sect.|3| The analysis of the limits 
IujIk. ^ in those sections is based on approximate calculations. In Sect. 0] we study 
the sticky limiting behavior of the trajectories obtained this way in more detail and we 
analyze its implications for the tracer transport. To confirm the approximate analysis of 
the first sections, we rigorously show in Sect. Elthat, in the intermediate compressibility 
regime, the stationary 2-point structure function of the forced tracer indeed converges 
when ~> in the fine-tuned way to the one corresponding to the sticky behavior 

of Lagrangian trajectories. Finally, in Conclusions, we summarize the obtained results 
and mention some open problems that they raise. Appendices contain more technical 
material relevant for the main text. 
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2 The stochastic process viewpoint 
2.1 Natural scale and speed measure 

In the intermediate compressibility regime, zero is a regular boundary for the unregu- 
larized operator M of Eq. H17() according to the Feller criteria. That is, if M is viewed 
as the generator of a stochastic process r{t) on then the realizations of the process 
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may hit zero and come back from it, with positive probabihty. In order to completely 
describe the process, one needs to specify a boundary condition for M at zero. On the 
other hand, for l^^lf^ finite, the regularized generator M^^^^ of Eq. (|13() should be taken 
with Neumann boundary condition at zero. As we have discussed above, this is because 
the problem on the half- line arose as the rotationally invariant sector of a non-degenerate 
problem defined on the d-dimensional space. 

To see what boundary condition is obtained when 0, it is useful to study the 

natural scale and the speed measure (see Appendix^ of the process r{t) with positive 
IvjIk and to determine their behavior in the aforementioned limit. The natural scale 
S{r) is the (positively oriented) coordinate in which the generator of the diffusion r{t) 
is without drift so that S{r{t)) is a martingale (there is of course an equivalence class 
of such coordinates, related by affine transformations). S{r) is strictly increasing and 
continuous in r, so that the function s{[ri,r2\) = S{r2) — S{ri) defined on intervals can 
be extended to a measure on absolutely continuous w.r.t. the Lebesgue measure). 
The measure ds is sometimes called the natural scale measure. Often it is more practical 
to work with the density of a measure than the measure itself. For a measure dfj, and a 
coordinate r on M_|_ we denote by fi{r) the density of d/i with respect to r. In the case 
of ds we have in particular s(r) = ^^§7^- 

In terms of the functions fi{r), f2{r) entering the generator M^^^ one can choose for 
the density with respect to r of the natural scale measure 

where ro > is an arbitrary point of the open half- line M^J,. It is easily checked that in 
the corresponding coordinate S, the process becomes driftless (no first order derivative): 

M,,^ = Doh{r)s{rfdl = Dom{S)-^dl. 

We shall define the speed measure as the measure dm with density 771(5*) with respect 
to the coordinate S. For convenience, we have multiplied the speed measure by 2Do 
with respect to the conventions of Appendix ^ in order to avoid such factors in later 
expressions. The diffusion process S{t) — S{r(t)) is then a (random, i.e. realization 
dependent) re-parametrization of the Brownian motion W{t), i.e. S{t) — W{T{t)), with 

2-Do ^ = m{S)\s=w(r) ■ 

For later use we also calculate the density of the speed measure with respect to the 
coordinate r: 

mir)^MS)f - ifiirMr)]-' = cxp ^£^^dr'^ (23) 

with the same tq as for the natural scale above. 

We expect that the process obtained in the l^Ji, limit corresponds to the 
natural scale and the speed measure that are appropriate limits of the same objects for 
Ik, Iv finite. In the next subsections, we shall analyze those limits in an approximate way. 

2.2 Approximate calculations 

We should know the dependence of functions /i and /2 on the regularization scales Iu,Ik,- 
If li, < Ik. then the smoothing of the velocity field is only significant at scales where it 
is already not the advection but the diffusion term that dominates, that is to say, the 
smoothing of the velocity plays no role and we may take 1^ = immediately. This case 
goes without difficulties and has been studied in ^U], see also [S]. The case which is 
interesting for us is the other one, when 1^ < lu- In this case we feel the diffusivity 
between and , then we feel the smoothing of the velocity field between Z„ and and 
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finally above l^, we are in the genuine Kraichnan regime with some scaling exponent ^ of 
the second order velocity structure function. 

In the present section, the calculations will be made by replacing functions /i , /2 by 
functions glued piecewise from pure powers, representing the different scaling behaviors in 
different subintervals. The gluing is done so that the function /i stays continuous. That 
the model functions obtained this way are indeed correct approximations of functions 
/i, /2 for the regularized Kraichnan model is shown in Appendix^ Thus for r > l^^ we 
shall take and /2(r) as in the (scale-invariant) Kraichnan model, that is — 

and f2{r) = a^r~^fi{r) with defined in (fH^ . For < r < the functions /i,/2 
will behave as in the smooth Kraichnan model (i.e. for ^ = 2), and we have to match 
them at li, with the values already given above (approximately; what matters, as we shall 
see, is only the order of magnitude and the ratio of the functions). Thus we shall take 
fi{r) — r'^ll^^ and /2('") = a2i"~^ fi{r). In the same way for r < the functions /i and 
/2 will be like in the pure diffusive case which, incidentally, is the same as the Kraichnan 
model with ^ = 0. With the same kind of matching as above, but now at Z^, we get 
/i(r) = and f2{'r) — {d — l)r^^/i(r) there. Table|21sums up our choices for the 

approximate versions of /i and /2. 



Scale 


Mr) 


rf2{r)/h{r) 


re [0,L] 




d-l 








r £ [liy, oo] 


r« 


a? 



Table 2: Approximate versions of functions /i and /2 

With the pure power choices for /i and /2, it is straightforward to calculate the 
natural scale and the speed measure on the intervals [0, ^k], [^k, In] and [Z^, oo). Let us first 
evaluate the ubiquitous subexpression s(r) = exp(— f2{r')/ fiir') dr'). We see that in 
each of the above intervals, f2(j'')/ fi{r') = a/r' with some a constant on the interval, so 
the integral evaluates to logarithms. Again, results are summed up in Table |21 



Scale 


s{r) 


m{r) 


r€[0,L] 












r G [li,, co] 







Table 3: Densities of the natural scale and speed measures 



2.3 Limit of the natural scale 

It is easy to see that if we maintain rg > fixed, then the density sir) of the natural 
scale measure tends pointwise to 

so{r) = rl'r~'^^ . (24) 

Recall that s{r) = ^^§7^- Thus S{r) has to be the integral of s(r), but we have the 
freedom to choose the constant of integration. We adopt the choice 

S{r)= / so(r')dr'+ / s(r')dr'. 

Jo Jro 

This definition has the advantage that, for r > l^, the value of S{r) is independent of 
Ifc^lu, because s(r) — so(r) there. In fact we have 

Sir) = So{r) = I so{r') dr' = —^r^'r^-'^i if r > . (25) 
Jo 1 - 
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In particular, in the limit I^tIv 0, the natural scale function becomes equal to So{r) 
for r > 0. Note that 50(0) — 0, observe, however, that if d > 2 then S{0) = —oo as 
long as Ik. is positive since in that case the integral s(r') dr' diverges as r goes to zero. 
Relation (|25|l may be inverted to give 

r(5) = [(1 - a^y-'^'S]^ if 5 > -^rlHl^'^' . (26) 



2.4 Limit of the speed measure 

Let us now calculate the speed measure of each of the intervals [0, Z^], [^k, I A and for the 
sake of completeness [l,^, R] for arbitrary R > l^. 

For [0, Z„] we get 




_ "[O ;(a5 + l-5) + (l-<»2),a2-l 

~ d " 



In the intermediate compressibility regime, + 1 — ^ > and 1 — 02 > 0, hence we have 
a positive power of li, and a negative power of l^- Depending on the way in which Zj,, 
go to zero, to([0, 1^]) can tend to zero, infinity or some fixed constant. For the finite limit 
one should have ;K+i-C)+(i-'>2) ^ ii-a^ ^ j 

In oc I, (27) 



The exponent of Z^ being greater than 1, this relation is compatible with the condition 
Zk ^ III in the limit Z^,Z„ -^0. In a moment we shall explain how the limit Z^,^^ ~^ 
taken with condition (|27|l corresponds to the sticky boundary condition. For now let 
us simply give some equivalent formulations of H27|l . From 114|l . it can be seen that 
2k = Dofi{0). Within our current modeling of /i, we have /i(0) — so that 

K^'-Doli-^ ll. (28) 



We shall use this equation to fix the relation between k and Z^ also in the exact approach. 
Recalling definition we see that the Prandtl number Pr = ^ is proportional to 
Q-v/lnY foi' SO that imposing relation l(?7| is equivalent to 

K (X li, ^ or Pr (x ly ^ . 

In particular k goes to zero while Pr goes to infinity when l^, goes to zero. 

The calculation of the speed measure of the interval [1^ lu] is performed similarly: 



02-1 

Now, since in the intermediate compressibility regime 02 — 1 < 0, in the limit 1^,1,^ ^ Q 
it will be that dominates Z"^"^ so that 



'm{[ln,lu\) 



1-02 



,(a5 + l-4) + (l-a2),a2- 
Iv Ik. 



We are in exactly the same situation as above and the same conclusions hold. 
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Finally the speed measure of the interval [li,, R] for some arbitrary i? > is 
»fl rR / ^ \ 



m{\l^,R]) — / m{r) dr — / r^^ ( — ) dr 
J I. J I. VaJ 



and it tends to a finite limit when 1^,1^ — > 0. 

At this point we may describe explicitly the speed measure dmo obtained as the weak 
limit of speed measures for positive l,^, 1^ when Zj^, — *■ in such a way that 

;(a« + l-0 + (l-a2);a.-l ^ ^ (29) 

for some A € [0, +00] (note that A has dimension (length)"'^^^ ^). On the open half-line 
of r > the density of dwo is 

r7io(r) =ro"''«r''«-«. (30) 
Besides, dmo has a mass at zero, given by 

mo({0}) = hm (m([0, L]) + m([Z., M)) = ^ + ^ " "J r""^ A . (31) 

It is convenient to characterize the "stickiness" of the boundary at r = by the quantity 
A = mo({0})/[r^^°£TOo(^)]r=o which we shall call the glue parameter. We obtain 

A^^^f^A (32) 
(1 - 02)0! 

This relation will be compared to the value produced by another approach in the next 
section, and finally to the exact result H62(l calculated using the precise forms of /i,/2 
instead of the approximate versions of Table El 

2.5 Convergence on the natural scale 

Here we shall show that the operations of changing coordinates from the original scale r 
to natural scale S and of taking the weak limit of the speed measure dm when 1^,1^ ^ Q 
commute. Observe that for d = 1 the natural scale coincides with the original one so we 
shall be preoccupied here only by the case d> 2. 

In the preceding subsection we computed the limit of m(r), working on the original 
scale r. Here we shall pass first to the natural scale S and then take the limit of m{S) 
and finally check that the result agrees with that of the preceding subsection. This is in 
order to exclude that we are in the kind of pathological case presented in Appendix El A 
priori such a situation could arise here because the notions of weak convergence on the 
original scale and on the natural scale do not coincide. 

The density of the speed measure w.r.t. the natural scale S is m(S') — 'm,{r)s^^(r). 
From (|26|l it is clear that for any > we can take l^, sufficiently small to have r(S) > li, 
and thus m(r) = mo(r) and s{r) = so('") so that m{S) = mo{r)sQ^ {r) — mo(S'o). Hence 
the trivial convergence of the density of the speed measure on the open half-line of positive 
S to mo (5*0) in agreement with the result of the preceding subsection. 

We still have to see what happens for S" < and to show that the limit mo(5'o) of the 
density of the speed measure has a Dirac delta term at So = 0. Let us evaluate for 5* < 
the asymptotic behavior of r{S) as IkJi^ 0. Observe that S{l,y) = J^" so(r')dr' > 0. 
So if < then r{S) < 1^. Now 

S{1.) = r So[r')dr' - r s[r')dr' 
Jo Ju 



-'0 I'l' „ '0 '■I' 



l~a^ " l-a2 



1~ 



I \ 1-02 
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when with l,^ < li, smce both 1 — and 1 — a2 are positive. This shows that 

for S" < fixed, we must have r{S) < asymptoticaUy. Let us exphcitly treat the case 
d > 2. The case cZ = 2 is only a bit different because of the logarithmic divergences. For 

r < 



S{r) = S{1,) 



2-d 



s{r')dr' 



2-d 



„2-d 



Suppose that r — 0(1^) (consistency will be checked) so that ~ o{r'^ 
that S{Ik,) 0, we infer that, for fixed S < 0, 



Recalling 



r{S)^ {2-d)r^"' C'^^'ll^'^+'^'S 



id^2)r-''^-S) 



We immediately check that r{S) — o(/k) so that our estimates are consistent, as promised. 
It can also be seen that this asymptotics for r(S) is uniformly valid on (— oo, S] for any 
S < 0. Now we may write 



m{S) = [m{r)s ^{r)]r=r(s) 

^ \d-2)-^''-^Kl' [IJhY-'^Hl-''' 



We infer that m{S) decays as a power of —S when S — > — oo, fast enough to be integrable. 
We also see that the coefficient in front of the power of — S" goes to zero when ^ Q 
with <l^. 

We still have to show that mo(S') has a Dirac delta contribution at zero. This can 
now be seen because 



m{S) dS 



m{r) Ar 



and we have shown that the right hand side has the limit given by (EH)- On the other 
hand we have just seen that for any S < Q 

rS 

m{S')dS' — > 

and also that S{l,y) > and S{l,y) 0. This permits to conclude that mo{So) has a 
Dirac delta contribution at 5*0 = with the coefficient mo({0}). An analogous reasoning 
for the case d = 2 gives the same result regarding the (weak) convergence of the measure 
dm on the natural scale. The limit measure coincides with the limiting speed measure 
in the r coordinate re-expressed in the limiting natural scale coordinate Sq. 



2.6 Conjecture 

It is legitimate to expect that the process r{t) converges in an appropriate sense to the 
one corresponding to the limit natural scale function 5*0 > given in (|25|l . and the speed 
measure with density 

moiSo) = T-""^"' [(1 - a^)^o]"'^^ + mo({0}) <5(5o) . 

Such a process is, in the natural scale coordinate 5*0, a reparametrization of the Brownian 
motion |M^(t)| instantaneously reflecting at zero: So{t) ~ \W{T{t))\, see Appendix IXI 
Away from zero, 2Do ^ — 'TJo(S'o)|so=|H'(T(t))| • At zero, the reparametrized process 
spends time proportional to mo({0}) times the local time at zero of |M^(t)| (although 
never an uninterrupted open interval of time) . This is the "sticky" or "slowly reflecting" 
boundary behavior PIEH] with the glue parameter A cx mo({0}), see HH2|I . The extreme 
cases A = and A = oo correspond, respectively, to instantaneous reflection and total 
absorption of the process at zero. 

This was our first argument for the sticky boundary behavior. 
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3 The linear operator viewpoint 



Our next argument is based on the study of the (generahzed) eigenfunctions of the 
operator M^^^ of Eq. (fT^ in the hmit 1^,1^. — > 0. If we take /i,/2 composed as above 
from pure powers, then we are able to calculate exactly the eigenfunctions and trace 
their limiting behavior. This argument is somewhat more shaky because it needs some 
explicit, not totally realistic (though representative), form of the functions /i and /2 
(before, such a representation was really required to hold only approximately). 



3.1 Derivation of the boundary condition at zero 

Let us consider the differential operator D^'^My ,^ — fi{r)d^ + f2{'r)dr with the approxi- 
mate form of functions /i and /2 given in Table|21 Thus for r > 1^ we take /i (r) = and 
/2(r) = a^r~^ fi{r) so on the interval [1^, oo] a pair of linearly independent eigenfunctions 
corresponding to the eigenvalue —E of D^^AIiy^^ is 

^Ur)^r'-^J±b,{VEr"-^) 

where 6^ is given by H21|l and Jb is the Bessel function of the first kind of parameter b. Note 
that these eigenfunctions are independent of 1,j,Ik. When we impose at zero a boundary 
condition^ in general the eigenfunctions of the operator with boundary condition will 
form a one dimensional linear subspace of the two dimensional linear space spanned by 
That is, if the linear combination c^$^ + '^e^e eigenfunction verifying the 

boundary condition at zero, then c^/c^ is fixed (depending on E and the boundary 
condition, of course). In part we shall proceed in the opposite way. First we calculate 
c^/c^ for Ii/^Ik finite. Next we deduce the effective boundary condition at zero that 
would give the same quotient if 1^,1^. were zero. Finally we look at the limit of c^/c^ 
when Zk, 0, and deduce the "limit" of the effective boundary condition. 

First of all we give explicitly a generating pair of eigenfunctions of D^^Mi^j^ asso- 
ciated with the eigenvalue ~E, on each of the intervals [0, [Ik,Iv], [liy,oo]. We shall 
denote these six functions gi , , respectively. Because at zero we have the Neumann 
boundary condition, on [0, 1^] we give only the corresponding eigenfunction: 

g:[{r) ^Clr^J_^i{^/EJ^r) with Ci = 2'^! (^1 - [E/ k)^ 

(the normalization was chosen for later convenience) . On [Z„ , we may set 

1 - as ± a/(1 - 02)2 - All^^E 
g±(r-)=r^±(^.'^) where 7±(i^,M = 2 

(recall that 1 — 02 > in the intermediate compressibility regime). Finally on [Zi/,00] we 
take 

53^(r) = 2±''«r(l±6e)i?T^<i>|(r). 

In order to construct the eigenfunction of M^^^. on the whole half- line [0, 00], we have 
to find the correct linear combinations of the generating pairs on each sub-interval by 
matching the functions and their derivatives at each border point. Formally, if we have 
to match c^g^ + c^.9j^ with cj'g^^ then, in matrix notation, we should have 

9t'{r) '{r)J J \gf{r) g^ ' (r) 
where g'{r) = drg{r). That is to say. 




'^Exactly one, given as LI^b] = where L is an operator that can be expressed as the limit at 
zero of some finite linear combination of derivatives of order or higher, with not necessarily constant 
coefficients . 
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Now, because on [0, 1^] we take as the eigenfunction, it follows that, by making the 
above described matchings, we get on [1,^, oo] the linear combination c^g^ + c^g^ with 



'E 



c^eJ Kgtiiu) 93 ' {^u)) \gt {lu) 92 (^•y)J Kgt'i^n) 92 ' (in) J V.9l '(^«) 

Note that 



(33) 



5+ _ T{l + b^) ( 4 



c. 



For a general matrix ( ° 5 ) its inverse is } u \ ''1 ■ We use this formula to rewrite 

yc dj ad— bey— c a J 

(|33|) . but since we are only interested in the ratio c^/c^, we get rid of the factor. 



Thus 



93 -93 m\ f9t{h) 92 C-f) 

-gfil.) 9t{lu) ) \9t\l.) 92'{Q 



92 '(^k) -92 (^k) \ f 9l (Ik) 



(34) 



This gives an explicit formula for c^/c^ for l^^li^ finite. In Appendix lUl we show that 
if we take li,, 1^ to zero in a way that the limit (|29|l exists (notation Zy, — * ) then 

r d+l-a2 , , 

Imi = — 71>'E. (35) 

L,i.^o Ce [l - a2){l - a()d 

It is not difficult to identify the boundary condition implied by this relation. Recall 
that the eigenfunctions of D'^^M^^^ associated to the eigenvalue —E are of the form 
^b('') = cuffs' (^) + c^(?;^(7') for r > ly, (that is for r > if IutIk 0). Furthermore we 
see that 

g3-(0) = l, r'^^g-\r)\r^^ = Q, 

53+(0) = 0, r'^'gt' ir)\r=o = l-H- 

From this we conclude that verifies 

$b(0)=c^ and r'^«$'s(r)|,=o- (1-05)4- (36) 

For convenience let us denote Dor'^^dr = F and, as before, A/q.o = M. Now H36|) with 
imply that 

Da\F<i>Em = ~~XE<i>EiO) 

where A is given by (|32|l . Since $£; is an eigenvector of Dq^M with eigenvalue —E, we 
obtain 

(F$b)(0) = A(Af$£)(0). (37) 

This last form may be read as an independent of E boundary condition that all ^e verify. 
It becomes then the boundary condition for M. Note that the coefficient A is in [0, +oo]. 
We introduce thus the operator which is the version of M with the domain 'D{M^) 
corresponding to the boundary condition (p?7| . In the next subsection we shall show that 
is the Kolmogorov backwards generator of a sticky process whose glue parameter is 
equal to A. 



3.2 Interpretation of the boundary condition 

To see that boundary condition l|177|) corresponds to a slowly reflecting boundary, we may 
calculate the mass at zero of the speed measure of the diffusion process generated by . 
The speed measure may be defined, up to a multiplicative constant, as the measure dm 
with respect to which the generator is self-adjoint 21 . It is straightforward to see 



15 



that for r > we may take for the density of the speed measure with respect to r (recall 
notation of Sect. ^ 

m{r) = r"«"^ . (38) 

We have in particular the relation m{r)M = Dodr o r"''^ dr ■ Here, as in the following, if we 
write M instead of M^, we mean not the operator with a specified domain but only the 
corresponding differential expression. Let us denote the scalar product with respect 
to the measure dm by (. , .),„. For to be symmetric with respect to dm, we need that 
{gi, Mg2)m — {Mgi,g2)m for any 51,52 G T^iM^). Using for the density of dm on 
the open half- line (0, cxd), and putting in a separate term the effect of a possible mass of 
dm at zero, this may be written as 

/■oo 

5i(0)(M52)(0)m({0})+ / g^{r){Mg2){r)m,{r)dr ^ 

(M5i)(0)52(0)m({0})+ / (M5i)(r-)52(r-) m(r) dr (39) 
Jo+ 

where the limit 0+ in the integrals indicates that the integral is over the open half-line 
(0, cxd). The integral on the left hand side may be rewritten using integration by parts as 

/"OO 

/ 5i(r)(M52)(r)m(r)dr = 
Jo+ 

51 (0)(i^52)(0)-(i^5i) (0)52(0)+ / (M5i)(r-)52(r)m(r)dr. 

We then see that (|39|l holds if and only if 

[(Af5i)(0)m({0}) - (F5i)(0)] 52(0) - 51 (0) [(Af52)(0) m({0}) - (F52)(0)] 

(strictly speaking the above expression is defined in the limit r ^ 0). This condition may 
only be verified for all 51,52 € V{M~^) if, for all 5 G V{M-^), ([m({0})M - F] g){0) = 0. 
Recalling the boundary condition (|37() . this amounts to the equality 

m({0}) = A, (40) 

We infer this way that the speed measure of Afj^ has a finite mass at zero, which means 
that Mj^ describes a process that is slowly reflecting at zero. From (|4(J|I and H38|l we see 
that the glue parameter of the slow reflection is A. 

More generally, the speed measure found in this section is coherent with the one we got 
from the stochastic process treatment. Compare H30II . (|31|l with H38|l . 14()|l . respectively, 
to see that the two speed measures are proportional (the second equals rg^ times the 
first one). 

4 The sticky process 

Now that we know the effective time evolution operator on large scales, we may inquire 
how the corresponding transition probabilities of the Lagrangian dispersion and other 
related quantities look like, what is their long time behavior and the induced effect on 
tracer transport. 

4.1 Transition probabilities 

The transition probabilities P|(r; dr') for the process with generator Mj^ are given by the 
self-adjoint exponential operators e**^^ defined in L^(m(dr)). It will be, however, more 
convenient to use the more standard L^{dr) conventions for the kernels and explicitly 
write the contribution due to the atomic term of dm. In such notation, 

P|(r;dr') = e**^A(r, r') dr' + A [e**^Mr, p) p?-'^^] S{r')dr'. (41) 
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The term concentrated at zero assures the conservation of probabihty / P|(r;dr') = 1 
as is shown in Appendix IeI Its presence imphes that sticky trajectories spend together 
positive time. In particular, the mean of the time that two trajectories starting at moment 
zero at distance r spend together up to moment t is A Jg [e'**^^ (r, p) p^^"'^]p=o ds. 

To compute the kernel e**^^ (r, r'), it will be more convenient to work with Schrodinger 
type operators. To this end, we shall first change variables by setting u = jZ^^cT^^^'"^ ■ 
Recall from H20|l that this transforms M to the generator of a Bessel process of parameter 
~b with b given by (|21|l (to ease the notation, we drop the subscript ^). After passing to 
the new variable, we conjugate M by the multiplication by u^~'^ and finally multiply it 
by —1 to obtain the differential operator N of the form 



The boundary condition 1371) imposed on Mj^ becomes for N the condition 

u^'^''duu'-^-^{u)\u^a = ~^lC-^u''-"-N^{^l)U=,, (42) 
where Ch = 2^^^^ r{i-b) introduced for convenience, and 

_ (2-g)2^-ir(b) , 
^ - r(i-6) ^0 ^ 

is a positive constant proportional to A. We shall denote by A'^ the operator N with this 
boundary condition. Note for later use the relationship 

e*'^''~^ir,r')dr' = uirf-^ e-''^'^iu{r),u{r'))u{r')^-^du{r') (43) 

The spectral properties of the operator A^^ are studied in detail in Appendix IdI Here 
we rely on those results. The spectrum of A'^ is R+. The (generalized) eigenfunction 
lPi^_e{u) of this operator, associated to the eigenvalue E >0, may be taken as 

Vt.Au) = \J-^,(^fEu) - AiPi-'' Jb(V^u)] . (44) 

The (scalar) spectral measure corresponding to this choice (of normalization) of the 
eigenfunctions is 

dE 

^""^^^^ " 2(1- 2/1^1-" cos(fe7r)+Ai2s2(i-b)) ■ (45) 

The kernel of the exponential of the operator N may be written using spectral calculus 
as 

/>oo 

exp{-tN^){u,v) = e'*^^^,E{u)^,,,E{v)diy,,{E). (46) 
Jo 

An important aspect of the flow is the behavior of two particles released at the same 
point at the same time. The evolution of the inter-particle distance of such a pair is 
described by the r ^ limit of P|(r;dr'). Substituting (|46(l into (|43|l and that in turn 
into l|lT|l . we obtain: 

e-*^ E-i V,A<r')) du,{E)) u{r')i-' du{r') 

(47) 

+ '''-Z:f"'' ~x{fe--^E-^dME))sir')dr'. 

The first term on the right hand side is the regular contribution absolutely continuous 
w.r.t. dr' . It describes the probability that a pair of particles starting together at time 
zero are separated at time t by some finite distance r' with dr' precision. Presence 
of such a term indicates that the sticky Lagrangian flow should be stochastic in each 
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velocity realization, just as the instantaneously reflecting flow corresponding to A = 0. 
The second term on the right hand side of H47|) is concentrated at r' = and describes the 
probability that two trajectories starting together will be together at time t. It is absent 
for A = 0. Recall, however, from JT, that in the strongly compressible phase p > ^ one 
has lim P~{r;dr') = S{r')dr' signaling that the Lagrangian flow is deterministic there. 

\ : : . ■ : ■ : of both regular and singular contributions to 1471) in the sticky flow is one 
of the indications of a hybrid nature of this case. 

In the special instance of /i = or /i = oo, corresponding to instantaneously reflecting 
or absorbing boundary at r = 0, the integral in (|46|l may be calculated explicitly. We 
make the change of variables E = and then use Weber's formula, Eq. (6.633.2) of |12| . 
to obtain 

exp{-tNo){u,v) = / e-^'^uij-h{VEu)vij^b{y^v) — 
Jo 2 

= u^v^ J e J_6(zu)J_6(zw)zdz= -^^exp ( — ]I^b[^) 

where /_t, is the modified Bessel function of the first kind of index ~b. Analogously, 

dE 



exp{-tNoo)iu,v) = I e~*'^u^Jb{^/Eu)v^Jb{^/Ev)'^ 







cxp — — h 



2t V 4t J \2t 
We then have 

exp(iMo.)(r, r') = r'"^ exp Tlfllllfli^ j^, ( 2^zr±_\ (4^) 

It can be easily shown that these results are in agreement with the well known transition 
probabilities of the Bessel process with reflecting or absorbing boundary at zero [2]. 

Some other quantities of interest are the hitting times Hr' at some point r', in partic- 
ular, the doubling and halving times of the inter-particle distance, see e.g. 0], Sect. 2C. 
Their expectations are given by the following formulae (with notations of Appendix IXI) 



E.(e-°^-') ^ , if r<r' 



and 



,2D- 



where (jj^-a is the solution of the differential equation (iV -I- q)4) = that verifies the 
boundary condition 142(1 at zero and t^-a is the solution that tends to zero at infinity. 
We refer here to Appendix^and Appendix^ in particular to formula ((63|l of the former 
and formulae H68|l and ((69|l of the latter. 



4.2 Long time asymptotics 

The sticky transition probability densities P|(r; dr') of the inter-particle distance deter- 
mine the behavior of the passive tracer 2-point function obtained in the corresponding 
limit /k, — > 0. The free decay of the tracer 2-point function is described directly by the 
transition probability P|(r;dr'). Indeed, for an initial tracer distribution with a homo- 
geneous isotropic 2-point correlation function i^(0, r), at time t the tracer distribution 
will be 

F{t,r)^ / F{{),r')Pl{r-dr'), (49) 
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as follows from Q . Thus the long time decay of the tracer 2-point function is determined 
by the large t asymptotics of e**^^(r, r'). Similarly, with isotropic forcing as in (|HJ|) and, 
for simplicity, no tracer at time zero, 

pt POO 

F{t,r) = ds x{r')Pi{r;dr'), (50) 
Jo Jo 

as follows from 111(1 . Thus the long time behavior of the forced tracer 2-point function 
is determined by the large t asymptotics of e**^^^ (r, r') ds. We shall assume fast decay 
for large r of both F{0, r) in the unforced case and of x(r) in the forced one. 

As before, it will be more convenient to study instead of the operator N^. To 
obtain the large t behavior of e~*^^'(u,w) or of Lf^{t;u,v) = /g e"'*^^' (w, u) ds, we may 
consider their Laplace transforms given for a > (where a is the variable conjugate to 
t) by the resolvent kernel 



/•oo 

Jo 



and by ^ (A^^ + a)~^{u, v), respectively. By the well known Tauberian-Abelian theorem, 
for any real numbers —1 < pn < ■■• < Pi < oo, the following behaviors of the function 
f{t) and of its Laplace transform f{a) are equivalent 

N 

(i) /(t) = c,tP' + o(iP~) near t = oo , 

1=1 

N 

(ii) /(a) = CiT{p, + l)a-P'-i + o(a-P"-i) near a = 0+ . 

i=l 



The resolvent (A'^, + a)~^ of A"^ is studied in Appendix IdI It can be written like in ((7()|l 
(except for the change of sign of a) and expanded in a for a small to all orders < 0. It 
is easy to see that, for /i < oo, in the numerator of H7U|) only the leading term is of order 
strictly smaller than zero and there is exactly one term of order zero. More exactly, one 
obtains: 

1 1 1 

— Cb {uv)2~''- — - - — min(?i, w) 3"'' max(u, w) 3+'' + o(a°) 

1 + fia^ " 

= Cb(ww)3-'' J2 (-A^)" a"^'"''^"'' - ^ min(u,w)3-Siiax(u,w)3+^ + o(a°). 

0<»<T^ 



This permits to infer the expansions 



E 



0<n<, 



r{b-n{l~b)) 



^fc-„(l-6)-l ^ o{t-^+') (51) 



for any e > 0, and 



0<n<^ 



b-n( 

i- min(u, w)^"'' max(u, w)^"'"^ + o(l) . (52) 



2b 



They give the large t asymptotics relevant for the study of the long time behavior of the 
tracer 2-point function in the decaying and forced regime, respectively. 



19 



4.3 Consequences for the tracer transport 

Substituting the expansion H51() into H43() and then into 141() . one infers from (|49() that 
the dominant terms in the free decay of the tracer 2-point function are r-independent: 

JTt^ — ('^-'S ^2fc-l nb-l^' ^b-n(l-b)-l 

F{t,r) - (— ) D, C, ^ r(6-n(l-5))* 

r'"'^ F(^^y)Ar' + AF(0,0)) + 0(^-1+^). 

Similarly, from (|52ll . (|43|l . (|41() and, finally, (|50|) . we infer that in the forced situation 
F(t,r) = 

00 



T{l + b-n{l-b)) 



- 3;;^ (r-^-'^ y r""-^ X{r') dr' + ^ r''-^x{r') dr' + A r^--^ x(0)) + o(l) . (53) 

It is interesting to notice that the number of terms appearing in this expansion is variable 
depending on the value of the integer part of j-^^. In particular if n(l — 6) — 6 = 0, i.e. 
6 = 1 - for some natural number n, then there is a supplementary term of order 
in the forced 2-point function. 

For A = 0; the result 1)53(1 reduces to the one worked out in 11 . The physics of the 
solutions with < A < 00 is not very different from the one for A = described in Sect. 
3.2 and 3.3 of The tracer "energy" with density ^6'^ is dissipated but, in the forced 
case, also building up in the constant mode growing like . Non-zero A brings subleading 
terms in this buildup proportional to ^''-"(1-*'). in the stationary state the mean tracer 
energy density (^0^) is infinite. The rate of pumping of the constant mode oc t^~^ goes, 
however, to zero. The tracer energy dissipation rate equal to —^MF{t)\r=Q approaches at 
long times the stationary value equal to the injection rate ^x(O) and the stationary state 
exhibits a direct energy cascade. The tracer 2-point structure function approaches 
the stationary form 

S2{t,r) = 2{F(t,Q) - F{t,r)) 



1 



r r 

0<r' <r" <r 



/'-«.'«e-?^(,')d/dr" + ^Ari-'^«X(0) • (54) 



For small r the integral on the right hand side scales as r^~^ while the term multiplying 
A is proportional to r^~°'^ . Since 1 — < 2 — ^, a non-zero A has an important effect. 
It changes the normal small r scaling cx r'^^^ of the 2-point structure function occurring 
for A = to the anomalous scaling proportional to the zero mode r^~°« of M. 

Let us also briefly mention the case A = 00, i.e. the one with absorbing boundary 
condition and coalescence of trajectories that was not covered above (the limits t 00 
and A ^ 00 do not commute!). Using 1)48(1 . the transition probabilities take in this case 
the form 

PUr;dr') = e**^-(r,r')dr' + [l - ^(6, ^-!±^) r(6)-i] <5(r') dr' , 

the same as in the strongly compressible phase, see Eq. (2.26) of jJJ. The coefRcient of 
S{r') is obtained simply by calculating the missing mass of e**^*^ (r, r'). For large time t, 
the 2-point function of the forced tracer becomes: 

Fit,r) = {t- (,_gp(,i,)r(i+.) ^o''t^"^^^'° + ,,-,n,-b, D^'r'-')x{0) 
+ -^—-^ ( / r' x(r') dr' + r'-^ / r' x{r') dr') + o{l) . 
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The behavior of this solution is similar to that of the strongly compressible phase analyzed 
in Sect. 3.4 of and quite different from the one for finite A. The scalar energy builds 
up linearly in time in the constant mode in an inverse cascade process and there is no 
persistent dissipation. At difference with the solution in the strongly compressible phase, 
however, the tracer 2-point structure function does not reach a stationary regime due to 
the contribution proportional to the zero mode r^~" growing in time like t^~^. 



5 Exact limit of the stationary 2-point tracer struc- 
ture function 

In the main part of the paper we have used in the calculations the approximate forms 
of functions /i and /2 instead of their exact versions. The question arises whether such 
calculations reproduce the actual behavior of the limiting dispersion process, up to finite 
renormalization of the glue parameter A. Besides, one would like to have exact expressions 
for some quantities, like the mass at zero of the limiting speed measure or the limiting 
stationary 2-point structure function of the tracer. It turns out that it is indeed possible 
to obtain precise formulae and that the asymptotic behavior of the integrals that we 
calculated approximately differs from the exact one by a (finite, non-zero) multiplicative 
renormalization. 

In what follows, we shall need detailed information about functions /i,2(?'; lu, Ik) and 
their ratio /2(?'; 'i/, ^k)//i(''; ^i^, ^k) that is described in Appendix IFl In particular we 
prove there the positivity for r > and the scaling properties 

/i(r;;,) = /«/i(^;l), f^{r;l,) ^ li'' hi{-; 1) (55) 

of functions fi.2{r;liy) = /i.2(?'; 'i/, 0). We also establish two decompositions. The first 
one: 

hip; 1) = Co p'f*ip) =Cop'+ p'ftip) (56) 

/2(p; 1) = 02 C-o pff{p) ^a2Cap + p'ftip) (57) 

with lim /*2(/o) — 1 'wiU be used for small p. The second one: 

hip; 1) = p« /<>(/,) ^ p« + p^-'f^{p) (58) 
hip; 1) = «c P^-'fiip) = P^-' + P^-'f?ip) (59) 

with lim f^ip) = 1 will be employed for large p. The coefficients 02 and ae are given 

by (|18l) . The above decompositions are used in Appendix |F] to establish bounds (|83|l 
to H85|l on the ratio /2(''; ^i^, ^k;)//i(''; 'i/, ^k)- Given such bounds, it is immediate to see 
that there exists a constant B' such that, for any r > and any < l,^ < < 1, the 
density s(r) of the natural scale measure given by H22(l is comprised between 1/B' times 
and B' times the approximate expressions listed in TableEl Using additionally estimates 
(|81|) to H82() from Appendix IfI one infers that the same statement holds for the density 
mir) = [fiir)sir)]~^ of the speed measure. It is now easy to show the convergence of 
s(r) and mir) to soir) and mgir) given by H24|l and H3()|l for r > 0. Controlling what 
happens around r — will be the main difficulty. 

After the above preparation, we pass to the main topic of this section: the proof of 
convergence, when l^^l^, — > 0, of the stationary 2-point structure function of the forced 
passively advected tracer. We shall show that, under the condition H29(l with A < 00, 
the limit exists and corresponds to the structure function computed directly a,t l^^Ji^ = 
with the sticky boundary condition at zero, see H54|l . Besides, we shall establish a precise 
relation between A and the glue parameter A of the boundary condition (|37|) . 

For positive, the stationary 2-point structure function is the unique solution 

vanishing at zero together with its first derivative of the equation Af^^„52 — 2x- It is 
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given by the relation: 



S2{r) ^ J- s{r")dr" i 
°° Jo Jo 



m{r') xi^') dr — -;— / s{r") m{r') xi"!"') dr dr' 

Do 



0<r'<r"<r 



The last integral may be cut up into six pieces according to the positions of r',r" with 
respect to Ik,Iu- The bounds given for s{r) and m(r) indicate that in each domain the 
dominant behavior of the contribution to 5*2 (r) is estimated correctly by using the rough 
forms in Table |21 up to multiplication by a constant bounded independently of 
and r. We get the following behaviors. The domain < r' < r" gives a contribution 
behaving at lowest order like r'^~^ for small r and like r^~°« for large r. Domains 
l^<r' < r" and r' <l^<l^ < r" both give terms of order Vi-"^? . 

Finally the contribution from domains r" < r' < 1^ is always subdominant compared to 
the former ones. This suggests that the limiting value of 5*2(7"; ly, 1^) will depend on the 
lt^~^ as Iu,Ik go to zero. The only task left is to calculate the limit 
A with A < 00. 



limit of I 



of 52(r; for f-^+^^-^Hl^-^ 

We already know that the domains with < r" < r' < li, will give vanishing con- 
tributions (supposing that the others give a finite one). Let us turn to the other three 
ones. The easiest to handle is ly < r" < r' < r. Because of the bounds on s(r) and m{r) 
we can immediately use the Dominated Convergence Theorem to obtain 



lim 

Iu,Ik^O 



s{r"; ly, Z„) m{r'; 1^, x(r') dr' dr" 

Ii/<.t' <ir" <ir 

s{r"- 0, 0) m{r'; 0, 0) xi^') dr' dr" = 

0<r' <r"<r 



r ^ ^ ^ ) dr dr . 

0<r' <r" <r 



Finally we turn to the two remaining domains covering the region < r' < < r" < r. 
The corresponding integral factorizes: 



/ 



s{r";ly,l^)m{r';ly,l^)x{r')dr'dr" 



0<r'<l^<r"<r 



s{r"-ly,l^)dr" 



m{r';lyj^)x{r')dr' 



To the first factor, one may again apply the Dominated Convergence Theorem to obtain 



lim 

l..U^OJi 



s{r"-ly,l^)dr" 



s(r";0,0)dr" = / ( ) ^ dr" = 



1 — 



Control of the second factor, that converges to the mass at zero of the limiting speed mea- 
sure multiplied by x(0), is the crucial element of the argument and it requires technical 
work. We postpone it to Appendix O We prove there that for A < oo, 



lim 

lu-l^-^O JO 



m{r';ly,h)x{r')dr' = r/«rAx(0) 



(60) 



with 



Y 



1 [l^o 



■ exp 



S(^,i^) 



/2*(p)-«2/r(p) 



Coff^P) 



dp 



dp 



(61) 
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where l^Jn^O denotes the hmit I^Jk, ^ under condition (|29|l and B{x,y) is the 
Euler Beta function. The constant Y is positive. The mass at zero of the hmiting 
speed measure is then equal to Tq '^Y X which should be contrasted with the approximate 
formula Given the above result, the contribution of the region < r' < l^, < r" < r 
to 52 (r) becomes 

2 Y 

Xr^-'"^ x(0). 

The final outcome is the relation: 

|-( J r"-'^^r"'^-«x(r')d/dr" + ^Ar^-^xW) 

0<r'<r"<r 

which agrees with expression H54|) of Sect. 14.31 if we put A = YX. This means that the 

structure function obtained in the Z„ ^ limit coincides with the one that is obtained 
directly using the sticky process with the generator Af^ for the value 

X = YX (62) 

of the glue parameter. Note that the effect of using the exact versions of /i, /2 instead 
of the approximate versions of table Table |21 manifests itself only in the change of the 
proportionality constant between A and the glue parameter A, i.e. in a finite multiplicative 
renormalization of A. 

In view of all that has been said, it is natural to expect that the entire Lagrangian 

pair dispersion process behaves in the 1^,Ik,-^Q limit as the diffusion process with the 
generator for the glue parameter A given by (|(j2|l . In particular, such a process is 
equivalent in the coordinate H19|l to the Bessel process of parameter ~b^,p with a sticky 
behavior at zero. 

6 Conclusions 

We have analyzed in this article the small viscosity, small diffusivity behavior of the 
Lagrangian dispersion in the Kraichnan model with intermediate compressibility degree 
+ ^ < p < In this interval, the Lagrangian trajectories may separate fast 
due to the spatial roughness of velocities but may also come close due to the trapping 
effects of compressibility. As first suggested in [5] E], we have discovered different possible 
asymptotic regimes of the Lagrangian flow, depending on the limiting behavior of the 
Prandtl number when the viscous and diffusive cut-off scales li, and are taken to 
zero. This arbitrariness reflects the frustration of the particles unable to choose between 
opposites trends of life. Specifically, we have argued that there are different limits of the 
dispersion process depending on the behavior of the combination 

of the cutoff scales. If this combination goes to zero when li, and are sent to zero, the 
resulting dispersion process is that of trajectories instantaneously reflecting off each other 
upon hitting. If it goes to infinity, the trajectories coalesce when they meet, behaving 
similarly as in the strongly compressible regime. Finally, if the above combination goes 
to a finite limit A when 1^,1^ — *■ (which sends to infinity the Prandtl number at a 
specific pace) then the resulting dispersion process exhibits the sticky or slowly reflecting 
behavior with A proportional to the amount of "glue" keeping the particles together. Such 
a behavior leaves a visible imprint on the passive advection of tracers in the subleading 
contributions to the tracer energy condensation in the zero wave number and, even more 
dramatically, by generating anomalous scaling of the stationary 2-point structure function 
of the tracer. The analysis in the main part of the paper was based on approximate 
calculations of the asymptotic behavior of the natural scale and speed measure of the 
Lagrangian dispersion process and of the eigenfunctions of its generator. To set the results 



lim S2{r) — 
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on a firmer ground wc have also shown rigorously that the finite l^, Z„ stationary 2-point 
structure function of the tracer converges, under the limit l,^, ^ with finite A, to the 
stationary structure function obtained directly from the sticky dispersion process. This 
argument allowed to fix exactly the proportionality constant between A and the value of 
the glue parameter. Undoubtedly, with a little more work controlling the convergence of 
resolvents of M^^^. to that of one should be able to prove, along the lines of 9 that 
the laws of the dispersion processes for positive 1^,1^ converge to the law of the sticky 
process. 

The main open problem, untouched by our analysis, is the construction of iV-particle 
processes corresponding to the sticky behavior of the two-particle dispersion. In partic- 
ular it would be interesting to know whether the amount of two-particle glue is the only 
parameter that labels possible Lagrangian flows in the moderately compressible phase 
of the Kraichnan model. The Dirichlet form approach used in |18j in the 1-dimensional 
^ = case to tackle such questions is unavailable in the other instances, at least in 
its classical form, due to the lack of symmetry of the generators of the A^-particle pro- 
cesses. Further open questions of fundamental importance concern possible occurrence 
of sticky Lagrangian flows in more realistic velocity ensembles exhibiting fully developed 
turbulence. 



Appendices 

A Briefly on one-dimensional diffusion processes 

We collect here some facts, used in the main text, about one-dimensional diffusion pro- 
cesses. The reader may wish to consult the relevant literature (e.g. [2ii2Q,) for an extensive 
treatment. 

Let X{t) be a stochastic diffusion process on R+. The Kolmogorov backwards evo- 
lution operator P* of the process, acting on bounded continuous functions / defined on 
R+, is given by 

{P\f){r) = Erf{X{t)) = I f{r')P\r,dr') 

where we define Erf{X{t)) = E{}{X{t)) \ X(0) = r). The measures P*(r,dr') giving 
the kernels of operators P* are the transition probabilities of the process. The family 
(P*)tgR^ forms a one-parameter semigroup whose generator dtP*\t=o we shall denote by 
M . In general M is a second order differential operator. The transition probabilities 
verify the PDE 

dtP\r, dr') = M{r)P\r, dr') . 

Any regular diffusion process may be derived from Brownian motion by an adequate 
change of variables defined in terms of the natural scale and the speed measure associated 
to the process. The natural scale is defined as the unique, up to affine transformations, 
continuous strictly increasing function S such that S{X{t)), i.e. X{t) considered in the 
new spatial coordinate S, is a martingale. S is called the natural scale and s(r) = ^^§7^ 
the density of the natural scale measure with respect to the coordinate r on M-|_ . 

If the range of S is the whole real line then the process S{X{t)), for X{0) fixed, 
has the same distribution as an appropriately time-changed Brownian motion starting 
from S{X{0)). This means that if is a Brownian motion starting from S'(X(0)), 
then S{X{t)) has the same law as W{T{t)) for some function T{t) — T{t;W) which 
depends also on the realization of the Brownian path. Let us denote by L{t, S; W) 
the local time of the Brownian path W at the point S up to instant r. Formally, 
P(t, S; W) = /J" S{S — W{a)) da. The speed measure of the process X{t) is defined as 
the unique positive measure dm{S) on R such that the relation between t and r may be 
written as 

T{t; W) = inf j cr : / L{a, S; W) dm{S) > t \ 
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or, equivalently, 

t{T;W) = / L{T,S;W)dm{S) 



Intuitively, the speed measure expresses how much time the (time-changed) Brownian 
motion needs to advance "one step" at a given point in space. The larger the (density 
of the) speed measure at some point, the slower the (time-changed) Brownian motion 
advances. 

The generator M of the Kolmogorov backward evolution semigroup may be written 
in terms of the densities m(S) — '^"^^^^ or m,[r) — '^"^''^ and s(r) as 

2m{S) 2m{r) s{r) 

Similarly, if the range of S* is a positive half- line, which may be chosen as M_|_, then 
analogous statements hold with the Brownian motion W replaced by the Brownian mo- 
tion reflecting at zero, i.e. by |VF(i)|. In any case, the behavior of the process at the 
boundary point ?' = may be classified according to Feller's criteria expressible in terms 
of the natural scale and the speed measure. When r = is a regular boundary point, 
the behavior of the process at this point depends on the its mass m({0}) w.r.t the speed 
measure and is reflected in the boundary condition for the generator M. 

Finally let us recall from jl3j the following formula for hitting probabilities of some 
process Xt. Denote by Hr' the hitting time of the process at some point r' . Then for 
a > 

E.(e-"«-) = |^^ (63) 

where (j) is the solution of (M — a)(j) = that verifies the correct boundary condition 
either on the left (i.e. at if the process is defined on the half-line) if r < r' or on the 
right (i.e. at infinity if the process is defined on the half- line) if r > r'. 



B A pathological case of convergence 

Let L be an arbitrary length and I a small scale that we shall send to zero. A diffusion 
process may be specified by giving its natural scale and its speed measure. Consider such 
a process on M+ whose natural scale is 



Sir) 



l-L+fr for re [0, /] , 



r for r G [I, oo) 

so that the density of the natural scale measure w.r.t. r is 



sir) = 



- J I 



^ for r £ [0, 1] , 
1 for r e [I, oo) 



Let us take for (the density w.r.t. r of) the speed measure 

, , Jf for re[OJ], 
m(r) = < 

I 1 for r G [I, oo) . 

We see that taking the limit I for r > we get so{r) — 1 corresponding to So{r) = r, 
however lim S'(O) = —L. This is not very different from the case studied in the Sect.|21 

where .^(O) — — oo as long as there is finite regularization. For the limit of the speed 
measure we get too('') = LSir) + 1. 

Let us now pass to the natural scale first. Then m(5) — = 1 for all S. Thus 
moiS) = 1. This is of course incompatible with the previous result. 
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The explanation of this phenomenon is the following. To the right of I the process 
is just Brownian motion. To the left of / it is Brownian motion on a segment of length 
L (with reflecting left end) "squeezed" into the segment [0,/]. Thus when we take I to 
zero, the limiting process will be a Brownian motion on ]R_|_ which "sees" an additional 
segment of length L to its left. This is of course a non-Mar kovian boundary condition. 

C Asymptotic behavior of eigenfunctions 

Starting from formula (|34|) we calculate the limit of c^/c^ when Ik,Ii, — > 0. We assume 
that Ik < lu, but not necessarily l^, <^ lu- The calculation can be done in the usual way 
by expanding c^/cj^ in a multivariate power series. First we expand every term in (|34|l 
to the minimal orders in Iu^Ik- That is, we expand the expressions into a power series 
(and possibly powers of logarithms) in l^, 1^ and keep all the terms such that there is no 
term of smaller order simultaneously in and in 1^. 

We list below the expansion to minimal orders in Z^, ^i/ of all terms appearing in 134() : 

1-02 

1 - a2 

9t{i.) = (i.V^ - i]""' gt'iL) = j+{Lr^-' ^ (1 - «2)C'^^ 

1-02 

g+{l,) = {l,y^ ^ ll-- g+\Q = 7+(/,)^^-l ~ (1 - a,)l-^^- 

<?3+(;.) ~ ff3+'(/.) ^ (1 - ae)/r« 

Note that (/«)''' cannot be further expanded without any additional hypothesis on the 
relative behaviors of l^, and 1^ because both the base and the exponent go to zero and 
one depends on 1^ the other on li,. All we can say is that 1% — o{{Ik)^ ) for any e > 0, 
and {Ik)'^ < 1 when 1^ < 1. We shall keep (1^)^ as it is in the expansions. 

Calculations to minimal order with expansions into multivariate power series are a 
little trickier then with univariate expansions. Even if the coefficients of some but not all 
monomials in the expansion simplify to zero at intermediate stages then possibly higher 
order terms should be taken into account because they could give rise to terms of minimal 
order^. This precaution is implicit in the computations. However, as long as there is 
only one term of minimal order (and it does not simplify out), no special care is needed. 

After having determined above the behavior of each term, we may now calculate the 
sub-products of (|34() to minimal order. First evaluate 

-93 (.^•^) \ ( dtiQ 92(1^) 
9til^) I \9t'iQ 92'iQ 

-(1-02)/-^^^ 




(2-0" 



Ja^-a2)^i "" -(l-a^)/^"'^ 



It is straightforward to check that every coefficient is different from zero, except for 
— 02 which may be equal to zero. So there the expansion has to be pushed further, but 

•^Take for example the sum of two polynomials {X + XY + Y^) + {—X-\-Y'^). The sum of the minimal 
order expansions is (X + Y^) + {—X + Y^) = 2Y-^ but the minimal order expansion of the sum is of 
course XY + 2Y^. 
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since each following term is of higher order both in l^, and in , we may simply replace 
ji-a^-a, by 0(/i-'^^-'^«). Next evaluate 



92 '('«^) -92 f 9i (L) 



-92 .92 7 \9i 



Again, it is straightforward to check that every coefRcient is different from zero. We may 
multiply together (exactly) the preceding two subexpressions to arrive at 



1 



{2-ifil-a2) 
4(l-?+a5) 



O ( ) + (1 - «2)(1 - a«) Z;"^/.-"^^ 
Now all the terms are of minimal order. Additionally, we have 



3— 02 — — 4 



(64) 



(65) 
(66) 



Let us first suppose that = o{ly"'^l~'''^). In this case c^/c~j^ goes 

to zero. This is so even if in H64I) the expansion for cancels out, meaning that the 
expansion should be pushed further. Indeed, subsequent terms in the expansion are of 
higher order thus asymptotically smaller. On the other hand, there can be no cancellation 
in the expansion of Cj^ because (|65ll and the hypothesis of the present paragraph combine 

Conversely, if we suppose Iv'^'^l^"'^ = o(Z^^'5^"^Z^^(Zk)'' ), then c^/c]^ goes to 00. 
This even if in H64|) the expansion for cancels out, meaning that the expansion should 
be pushed further. Once again, subsequent terms in the development are of higher order 
thus asymptotically smaller. And now it is in the development of that there can be 
no cancellation because (I66II and the hypothesis of the present paragraph combine to 

In consequence, if c\lc^ should have a finite non-zero limit, then I'^'^^^'^^l'^^il^)'^ 
and lu '^^ must be of the same order. In particular they are much bigger than the other 
two terms so there are no cancellations in ()64|l . Also it is easy to see that in this case 
cannot decrease faster than some power of l^- If we write 1^ — 0{11^^^"'^1~^{Ik)^ ) 

2-S-a2 + °; 

then, using (1^)'^' = 0(1), we have = 0{l^). This implies {l^)^' - 1. We 

may then conclude that 

d + 1- a2 p .2-^-02+05,02-1 
(l-a2)(l-a5)d " ' 

This expression has a finite limit if and only if ll ^ "2+05^02-1 g^gg some finite limit 
A, just as in H29|) . We obtain this way formula H35|) . In fact our proof shows that H35fl is 
valid also if A is zero or infinite. 



D Spectral analysis of 

The (operator valued) spectral measure for the operator may be evaluated with the 
help of the formula 

where G^_q is the resolvent 
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and S is a (Borel) subset of R. The kernel G^^a{u, v) of G^.q is determined by demanding 
that for each v G K+, as a fmiction of u G R+, it is in the domain of and that it 
satisfies the equation 

{N{u)-a)G^,^^{u,v) =d{u~v). (67) 

To calculate Gp,q(m, v) we note that for u v we simply have {N{u) — Q)G^^a{u, v) = 0. 
Viewed as a function of u, G^_q(m, u) verifies at zero the boundary condition H42|l and 
should go to zero at infinity. A solution of {N^ — a)(t){u) = satisfying the correct 
boundary condition at zero is 



(68) 



Here I±i, is the modified Bessel function of parameter ±6. The square-root is taken with 
its principal definition. Similarly, the solution of {N — a)ijj{u) — decaying at infinity 
is 

■iPa{u) ^u^/^Kb{^/^u) , a^M+ (69) 

where Kj, is the modified Bessel function of the second kind of parameter b. With the 
correct matching at u = w to assure that (|67|l is satisfied, we obtain 



0f.,Q(")V'a(f) 



if u < V . 
if u > V 



(70) 



where w^^a is the Wronskian of and ipa, i-e. 

Wp.a = (l>'^,,a{z)iJa{z) - (/)^,a(z)'0^(z) 

which is independent of z. To evaluate the Wronskian, we may use the asymptotic 
expansions for z — > +oo 



/.(z) = ^(l + 0(i)). 



see Eq. 8.451.6,7 of |12,. The resuh is 



Wu 



Next we have to calculate the discontinuity in G^^a{u, v) along the cut for a € IR+. It 
will be convenient to write — (a ± zO) = e^'^^a. This will automatically give the correct 
determination of every function. For u < v 



G 



{u,v) 



1 



,1/2 



1 + /iai~''e=F'r«(i-b) 



Now we employ the general formula 



7r /_^(z) - /^(z) 

^i^[z) = 77 ■■ 

2 sm i/TT 



for V i'L , 



see Eq. 8.485 of [HI, to get 



1 



+ /ie=F"(i-'')ai-''/b(e=F¥^u 
The relations 



-6)' 



,1/2 



,1/2. 



2 sin feTT 



I^e^z) = e^'^J,(z) , L{e-^z) = e-¥-J^(z) 
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following from Eqs.,8.406.1 and 8.476.1 of permit to move the determination from 
the function argument to a multiplicative coefficient. The final expression is 



1 



2 sin 67r 1 + uai-''e=F'ri(i-&) 



''JbiVau) 



It is now easy to calculate the discontinuity G^.a+io ~ Gf^^a-io and to get for a = E Cz 

^ 2(1 - 2^E^-^ cos6vr + ^^i^^d^.)) ^-^(") ^-^(«) (^1) 
where ip^^E{u) given by (|44|l is the (generalized) eigenfunction oi associated with the 



E Conservation of probability 

We shall show here that the transition measures P|(r;dr') of the sticky process given 
by H41|) are normalized, i.e. that / dr') = 1. This is clearly the case for i = since 
P-°(r; dr') = 5{r - r') dr'. But 



/ e'^'~^{r,r')dr' = / M\r')e'^'~^{r,r')dr' - -F{r') [e'^'^^{ry)m{r')-^]^,^^ 



dt 

where — m{r)M o m{r)^^ = dr o F o m{r)^^ with m{r) — r°'i^i and F ~ D^r'^^dr 
is the formal adjoint of M w.r.t. the L^{dr) scalar product. The expression on the right 
hand side has the interpretation of the mean fiux of the probability current through 
r' — 0. It should be balanced by the rate of change of probability to stay at r' = 0. 
Using the relations 

e**^^ (r, r') m(r')"^ = e*^^^ (/, r) m{r)-^ 

and the fact that e**^>'(r, r') satisfies as a function of r the boundary condition . we 
infer that 



d_ 
dt 



e'^^~>-{r,r')dr' = -~X M{r') e'^^'>- {r' ,r) m{r)-^ 

= -AAf(r')(e**'^M^, r')m(r')"^ 
= -Am(r')"^A/^(/)e**^H^ r') 



r'=0 
) r'=0 
r'=0 



= ^~X^m(rT'e'^'Hr, r') 
dt ^ ' ^ ' ^ 

It follows then from (|41|l that the time derivative of J P|(r;dr') vanishes so that the 
normalization of P|(r; dr') does not change in time. 



F Estimates of functions fi and /2 

We present bounds showing in which sense the approximations given for functions fi and 
/2//1 in Table 13 are correct. 

The functions /i,2(r; Zy, /«) were defined by ffT^ and l(T3|l with k given by Let 
us examine their behavior for k — setting fi.2{r]lu) = /i,2(r; 0). Let us first show 
that /i,2(r; lu) are positive for r > 0. Note that JH)) implies that 

r 2sin^ ^=:^ dk f ^ ^ 

d.,{r;Q ^ DoJ ^^j^ P(k,p)fmk\)^ ^ DoJ T,,{k)dk 
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where T{k) is a positive (r-dependent) matrix. On the other hand, (|14|l and (|15|l may 
be written as 

Mr) = I T,,(fc)^dfc, Mr) = r-'l T,,(fc) (5,, - ^) dfc 

for any vector r such that |r| — r. Positivity of T imphes 

(k) nr, > , (k) (% - ^) > 

and for r > it is not hard to see that strict equahties hold except for a set of k of 
measure zero. Hence positivity of fi,2{'i",lu)- 

In order to study the behavior of fi,2{r, Iv) for small and large r it will be convenient 
to express them in terms of the function 

/I _ ik-r ^ J? 
fn\k\)——. 

Note that g is a smooth even function (i.e. all its derivatives of odd order vanish at 0) of 
r > 0. By scaling, 

g{r;lv) = lig{^;l). 

Around zero 



g{p- 1) = C7op2 + 0(/) ^ith Co = / 4^ -P- 

Rewriting 



(72) 



we infer that around infinity 

g{p; 1) = CooP« + 0(p«-2) (73) 



with 



C'oo — 



1 _ e*p/p dfc r(^) 



The spatial covariance (i(r; ^,y) given by © may be expressed in terms of the function g: 



pd-l 
d-1 



[5,, - d^) r-'^ 1 (r') '^-i5(r'; L) dr'] . (74) 



Indeed, by rotational covariance, dij has to be a combination of Sij and To show 

that the coefficient functions are correctly represented above, it is enough to note that 
da and didij are, in virtue of ((HJ, expressible in terms of g and of its r-derivative g' as 
da — Dog and didij = Dop^g' and that the same relations may be recovered from H74fl . 
Now using definitions (|14() and H15() we obtain: 

^ l^) - ^ 9{r; I.) + ^ (gir; Q ' ^ [ r' g{r' ,1^) dr') , (75) 

%^j2{r-M - ^.9(r;/.) + f r'"-' g{r' -l,) dr' . (76) 

Functions fi.2{r\ly) inherit from g{r\l^) the scaling property H55|l . Relations H75() . (|76|l 
and expansions 1)72(1 and ((73|l permit now to write the decompositions 1(56(1 to ((59(1 with 
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Cq — ^^^CqDqI Dq and similar coefficient in (|58|) and (|59|l fixed to 1 by setting Dq = 

^d+i ^0 ■ Moreover the positivity of functions fi,2i'r] 1) and tlieir smoothness imply 
that there exists some constant B > such that for all < p < 1 



and for all 1 < p < oo, 



B 



< frip) < B 

B < ff{p) < B 



i = l,2, 
i = 1, 2 



(77) 
(78) 



B 



< fHp) < B 



-B<ff{p)<B 



i = l,2, 
i = 1,2. 



(79) 
(80) 



Let us come back to functions fi,2{f") lu, Ik) with < ^k, < 1- To obtain bounds for 
/i, we write for 1^ < r 



and infer that 



r^-<hir;L,L)<rHB + 2). 



(81) 



Similarly, the decomposition 

/i(r; L) = Co It'r'ff'i^) + 211-^1 

t'lf 

gives for 1^, < r < li, the bounds 



and for r <li^^ 



ll-^r^^ < h{r-l,,h) <li-\\C^B + 2) 



2li-X<fi{r;L,Q<li-^ll{CoB + 2). 



(82) 



This shows that /i(7';/^,/k) is bounded above and below by its approximate version of 
Table 12 multiplied or divided by a constant. 

Function /2 may be estimated similarly resulting in 0(1) bounds for r^ij j\. We 
need however (in particular in AppendixE} the following more precise estimates for that 
ratio: 



for < r : 

for Ik, < r < : 

for r < Ik. ■ 



f2{r\h,lK) 


_ 2± 




fi{r;li,,lK) 


r 




f2{r;L,lK) 


0.2 






r 




f2{r;L,lK) 


d- 


- 1 


fi{r;li,,lK) 


r 



7? 

< ^V^, 



< c 



I?, 



r 



(83) 
(84) 
(85) 



for some constant C which may be chosen independent of r, Ii^,Ik.- To establish H83|l we 
use the decompositions (|58f) and H59|) to write for < r 



f2{,r;l^,lK) a? 



rV2(f ; 1) + 2^li-Hl - ^ lihif-- 1) + 2li-Hl 



lihi^;l) + 2lt^ll 



f?if-)-^d?iT-) 



< 



B[l + \a^\] + 2\d-l-a^\ 1$ 



i_ 

B 
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where the inequaHty follows from the bounds (|79|) and (|80|) . In the same manner for 
< r < l^, using the decompositions (|56|) . H57|l and the bounds H77|) and (|78|l . we obtain 



Finally for r <l^, 



/2*(f)-«2/r(f 



^ r- I' K 



< 



B[l + \a2\] + 2\d - 1 - a2\ 

Co 



< 



r (7o|d~l-a2| + g[l + M-l|] 
Z2 2 



Let us end this Appendix by listing the scaling forms of the tensor d and of functions 
/i,2- Since g{r;0) = CocT^, we infer from (fTljl that 



dij{r,Q) = 



Dn 



d-1 



(i 



d + S. , ^{pd - 1) nrj 



and that fi{r; 0) = and f2{r; 0) — a^r^ ^, in agreement with the scaling form (|17|l of 
the generator M. 



G Exact speed measure at zero 

Here we prove the convergence (jHOJ, establishing the exact value of the mass at zero of 
the limiting speed measure. Using 123|l ). we have 



m{r';l^,lK)x{r')dr' = 



X{r') 



hir'^lu.h) 



exp 



fiir";l.,L) 



dr" I dr' . 



The integral in the exponential may be split into four terms: 



f2{r";L,L) 



'r' .fi{r";l.,L) 



dr" 



dr" 



.f2{r";l.,L) 



dr" 



We begin by calculating the last three terms as they do not depend on r'. The last one 
gives 



f2{r"X,l.) 



'"0 



dp 



and we know from (|83(l that the integrand on the right is dominated by Cp (^+5) (as 
soon as < 1), so we may use the Dominated Convergence Theorem to conclude that 



lim 



°° f2ip)-ae.f?{p) 



P'f^ip) 



dp. 
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As for the other two, 



a2 In ■ 



2 d—l — a2 
P 



Cop'ff'CtP) + 2 



dp 



with the integrand on the right dominated by C[(p)^p + p < 2Cp see (|51)l. so we 
may conclude that 



hm 



h{r";h,l.) 



dr" — 02 In \/ 7^ 



2(d - 1 - 02) , d - 1 - a2 , Co + 2 
■ dp — In ■ 



p(Cop2 + 2) 



Co 



Similarly, 



dr" 



a2 In 



/2*(p)-«2/f(p)]+2^^(f^) 



(7op2/i*(p) + 2(|f)' 



■dp 



where the integrand on the right is bounded by C[p + ij^)'^p < 2Cp, see (|84|) . Hence 



lim 



dr — 02 In 1/ — 



What remains to be evaluated is 

X{r') ( f2{r";l.,Q 



/2*(p)-«2/f(p) 



Coffip) 



dp . 



fiir';li.,L) 



exp 



xiLp) 



dr" dr' 



CoP^ff'Ctp) + 2 



exp 



i Coa2p'ftCtP') + ^^ ^ , 
CoP'^ff'CfP') + 2 ^ 



dp. 



Let us show that the Dominated Convergence Theorem applies once more. The following 
estimates are sufficient. For < p < 1, 

C'oa2p/2*(i:P) + ^ d-1 



x{Lp) 



c'op2/i*(bp) + i 



< X(0) , 



Cop'ftCtP) + 1 



< Cp. 



The first inequality is a consequence of xC*") < for r > 0, and of the positivity of 
/*. The second one is obtained by rewriting 1|85|1 using (j5t)|l and (|57|l . For 1 <C p < li/lln-i 
we use 



xiLp) 



Cop2/i*(hp) + l 



< 



X(0) 

Co 2 
B P 



Coa2pf*CfP) + 



2 



C'op2/*(j-p) + l 



02 



< c 



1 



Here, the first inequality follows from H77(l and the second one is obtained by rewriting 
(IHH) with the use of and lO. We then have 



exp 



.f2{r"-J,,Q 



dr" dr' 



x(o) 



Cop2 + 2 



exp 



iC'oa2p' + 2^ ; 

^dp 



Cop' 



/2 , 



1 



dp 



2^ (Co + 2)^^^ C-^ B{± i^) x(0) 
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where B{x, y) is the Euler Beta function. Gathering aU terms, in particular the powers of 
li, and Ik, that combine to ll""'^^ 12)^02-1^ obtain the result H6Q|I with Y given 

by 
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